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Abstract 

We present a novel approach to morph between two isometric poses of the same non-rigid object 
given as triangular meshes. We model the morphs as linear interpolations in a suitable shape space 
S. For triangulated 3D polygons, we prove that interpolating linearly in this shape space corresponds 
to the most isometric morph in . We then extend this shape space to arbitrary triangulations in 31^ 
using a heuristic approach and show the practical use of the approach using experiments. Furthermore, 
we discuss a modified shape space that is useful for isometric skeleton morphing. All of the newly 
presented approaches solve the morphing problem without the need to solve a minimization problem. 

1 Introduction 

Given two isometric poses of the same non-rigid object as triangular meshes S^^"^ and S^^"^ with known 
point-to-point correspondences, we aim to find a smooth isometric deformation between the poses. Inter- 
polating smoothly between two given poses is called morphing. We achieve this by finding shortest paths 
in an appropriate shape space similar to the approach by Kilian et al. 0|. We propose a novel shape space. 

A deformation of a shape represented by a triangular mesh is isometric if and only if all edge lengths 
are preserved during the deformation |9|. This property holds because each face of the mesh is a triangle. 
A deformation of a shape is called most isometric if the sum of the squared differences between the corre- 
sponding edge lengths of the two shapes is minimized. In this paper, we examine isometric deformations 
of general triangular meshes in 3D and of triangulated 3D polygons, which are triangular meshes with 
no interior vertices. We introduce a new shape space S for triangulated 3D polygons that has the property 
that interpolating linearly in shape space corresponds to the most isometric morph in R^. We then extend 
this shape space to arbitrary triangulations in 3D using a heuristic approach. Furthermore, we discuss a 
modification of the shape space that is useful for isometric skeleton morphing. 

2 Related Work 

Computing a smooth morph from one pose of a shape in two or three dimensions to another pose of the 
same shape has numerous applications. For example in computer graphics and computer animation this 
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problem has received considerable attention |[II1I3. A recent survey on this topic was written by Alexa Q. 
We only review the work most relevant to this paper. 

Before considering morphing three-dimensional mesh models, the two-dimensional version of the prob- 
lem has received considerable attention. In the case where the input is sampled over a regular domain, this 
problem is called image morphing. Image-morphing is widely studied and surveyed by Wolberg II2TII . In 
the case where the input is sampled over an irregular domain, the problem is to interpolate between two 
simple polygons in the plane. We only review work that makes use of intrinsic representations of the poly- 
gons. Sederberg et al. |[T4ll propose to interpolate an intrinsic representation of two-dimensional polygons, 
namely the edge lengths and interior angles of the polygon. Surazhsky and Gotsman llT9ll morph by com- 
puting mean value barycentric coordinates based on an intrinsic representation of triangulated polygons. 
This method is guaranteed to be intersection free. Iben et al. 1 8 1 morph planar polygons while guaranteeing 
that no self-intersections occur using an approach based on energy minimization. This approach can be 
constrained to be as isometric as possible. 

Sun et al. ifTSl morph between three-dimensional manifold meshes. They extend the approach by Seder- 
berg et al. fT4\ to three dimensions by extending the intrinsic representation to polyhedra. However, the 
developed methods are computationally expensive [IJ. We propose a different intrinsic representation of 
triangular meshes that offers the advantage of producing the most isometric morph and of being efficient. 

Sorkine and Alexa fTfj propose an algorithm to deform a surface based on a given triangular surface 
and updated positions of few feature points. The surface is modeled as a covering by overlapping cells. The 
deformation aims to deform each cell as rigidly as possible. The overlap is necessary to avoid stretching 
along cell boundaries. The deformation is based on minimizing a global non-linear energy function that 
is simple to implement. The energy is guaranteed to converge. However, since the energy function may 
have multiple minima, the algorithm is not guaranteed to find the global minimum. The approach tends to 
preserve the edge lengths of the triangular mesh. This property depends upon finding a global minimum of 
the energy function. One cannot guarantee to find this global minimum. 

Recently, Zhou et al. (24] proposed a new method to deform triangular meshes based on the Laplacian 
of a graph representing the volume of the triangular mesh. The method is shown to prevent volumetric 
details to change unnaturally. 

Recently, Kilian et al. |[9l used shape space representations to guide morphs and other more general 
deformations between shapes represented as triangular meshes. Each shape is represented by a point in a 
high-dimensional shape space and deformations are modeled as geodesies in shape space. The geodesic 
paths in shape space are found using an energy-minimization approach. Before Kilian et al. |9| presented 
the use of a shape space for shape deformation and exploration of triangular meshes, shape space rep- 
resentations were developed to deform shapes in different representations. Cheng et al. ||5l proposed an 
approach that deforms shapes given in skin representation, which is a union of spheres that are connected 
via blending patches of hyperboloids, with the help of a suitable shape space. Furthermore, algorithms for 
deforming curves with the help of shape space representations were proposed by Younes |[23l and Klassen 
et al. |[TOl . Eckstein et al. |6| propose a generalized gradient descent method similar to the approach by 
Kilian et al. that can be applied to deform triangular meshes. All of these approaches depend on solv- 
ing a highly non-linear optimization problem with many unknown variables using numerical solvers. It is 
therefore not guaranteed that the globally optimal solution is found. 

In this paper, we propose a novel shape space with the property that interpolating linearly in shape 
space approximates the most isometric morph in R^. For triangulated 3D polygons, we prove that the 
linear interpolation in shape space corresponds exactly to the most isometric morph in M^. For arbitrary 
triangulated manifolds in 3D, we provide a heuristic approach to find the morph. This heuristics is an 
extension of the approach developed for 3D polygons. The proposed methods do not require solving a 
minimization problem. 
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3 Theory of Shape Space for Triangulated 3D Polygons 



This section introduces a novel shape space for triangulated 3D polygons with the property that interpolat- 
ing linearly in shape space corresponds to the most isometric morph in R^. The dimensionality of the shape 
space is linear in the number of vertices of the deformed polygon. 

We start with two triangulated 3D polygons P^^^ and P^^^ corresponding to two near-isometric poses 
of the same non-rigid object. We assume that the point-to-point correspondence of the vertices P^^^ and 

p(i) 

are known. Furthermore, we assume that both p(^^ and P^^^ share the same underlying mesh structure 
M. Hence, we know the mesh structure M with two sets of ordered vertex coordinates V^^"^ and V^^"^ in 
R^, where M is an outer-planar graph. We will show that we can represent P^^^ and P^^^ as points p^^^ 
and p(^) in a shape space 5, such that each point p^^^ that is a linear interpolation between ^(^^ and p^^"^ 
corresponds to a triangular mesh P^^^ isometric to P(°) and P(i) in R^. 

As we know the point-to-point correspondence of the vertices P^^^ and P^^\ we can find the best rigid 
alignment of the two shapes by solving an overdetermined linear system of equations and by modifying 
the solution to ensure a valid rotation matrix using an approach similar to the one used for camera calibra- 
tion |[20l . 

Let M consist of n vertices. As M is a triangulation of a 3D polygon with n vertices, M has 2n — 3 
edges and n — 2 triangles. We assign an arbitrary but fixed order on the vertices, edges, and faces of M. The 
shape space S is defined as follows. The first 3 coordinates of a point p ^ S correspond to the coordinates of 
the first vertex v in M. Coordinates 4 and 5 of p correspond to the direction of the first edge of M incident 
to V in spherical coordinates. The next 2n — 3 coordinates of p are the lengths of the edges in M in order. 
The final 2(n — 2) coordinates of p describe the outer normal directions of the triangles in M in spherical 
coordinates, in order. Hence, the shape space S has dimension 5 + 2n — 3 + 2(n — 2) = 4n — 2 = 0(n). 

In the following, we prove that interpolating linearly between p(^) and P^^^ in shape space yields 
the most isometric morph. To interpolate linearly in shape space, we interpolate the edge lengths by a 
simple linear interpolation. That is, p^^^ = tp^j^"^ + (1 ~ ^ where p^^^ is the fcth coordinate of p^^\ 
The normal vectors are interpolated using geometric spherical linear interpolation (SLERP) |15|. That is, 
'P'k ^ ^^^sin'e^^ ^fe'^ + ^liae 'Pk^ ' whcrc 6 is the angle between the two directions that are interpolated. 

Note that the relative rigid alignment of P^^^ and P^^^ in R^ has an influence on the linear interpolation. 
That is, the interpolating shape space point varies as the relative rigid alignment of P^^^ and P^^^ in R^ 
changes. The change occurs because the angles between the normal vectors of p^'^^ change as a result of the 
rigid transformation. This is the reason we choose to find the best rigid alignment of P^^^ and P^^^ before 
transforming the polygons into S. 

To study interpolation in shape space, we make use of the dual graph D{M) of M. The dual graph 
D{M) has a node for each triangle of M. We denote the dual node corresponding to face / of M by D{f). 
Two nodes of D{M) are joined by an arc if the two corresponding triangles in M share an edge. We denote 
the dual arc corresponding to an edge e of M by D{e). Note that because M meshes a 3D polygon, it is an 
outer-planar triangular graph and so the dual graph of M is a binary tree. An example of a mesh M with 
its dual graph D(M) is shown in Figure [1} 

Theorem 1 Let M be the underlying mesh structure of the triangulated 3D polygons P^^^ and P^^\ The 
linear interpolation p^^^ between p^^^ andp^^^ in shape space S for < t < 1 has the following properties: 

1. The mesh P^^^ G R^ that corresponds to p^^^ G S is uniquely defined and has the underlying mesh 
structure M. We can compute this mesh using a traversal of the binary tree D(M) in 0(n) time. 

2. If P(^) and P^^^ are isometric, then P^^^ is isometric to and P(i). and are not 
perfectly isometric, then each edge length of P^^^ linearly interpolates between the corresponding 
edge lengths ofP'^^^ and P(^). 
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Figure 1: A mesh M with its dual graph D{M). 
3. The coordinates of the vertices ofP^^^ are a continuous function oft. 

Proof. Part 1: To prove uniqueness, we start by noting that the first vertex v of P^*) is uniquely determined 
by the first three coordinates of p^'^\ The direction of the first edge e of M incident to v is uniquely 
determined by coordinates 4 and 5 of p^^\ because each point on the unit sphere determines a unique 
direction in M^. The length of each edge of P^^^ is uniquely determined by the following 2n — 3 coordinates. 
Furthermore, the outer normal of each triangle is uniquely determined by the following 2(n— 2) coordinates, 
because each point on the unit sphere determines a unique direction in M^. Hence, the edge e is uniquely 
determined. For a triangle / containing e, we now know the position of two vertices of /, the plane 
containing /, and the three lengths of the edges of /. Assuming that the normal vectors in shape space 
represent right-hand rule counterclockwise traversals of each triangle, this uniquely determines the position 
of the last vertex of /. We can now determine the coordinate of each vertex of P^^^ uniquely by traversing 
D{M). We start the traversal of D{M) at D{f). Recall that the coordinates of the vertices of triangle / 
are known. Hence, when traversing an arc D{e) incident to we know the vertex coordinates of the 

shared edge between the two triangles corresponding to endpoints of D{e). Denote the endpoint of D{e) 
not corresponding to / by /^ For f , we now know the position of two vertices of f , the plane containing 
f , and the three lengths of the edges of /^ Hence, we can compute the position of the last vertex of f . 
Because we now know the coordinates of all the vertices of we can traverse all of the arcs in D{M) 
incident to D{f'). In this fashion, we can set all of the vertex coordinates of P^^^ by traversing D{M). 
Because D{M) is a tree, it is cycle-free. Hence, the coordinates of each vertex of P^^^ are set exactly once. 
Because the complexity of D{M) is 0(n), the algorithm terminates after 0{n) steps. 

It remains to prove that P^*) is a valid triangular mesh, that is, that the three edge lengths of each 
triangle of P^^^ satisfy the triangle inequality. We assume that both input meshes were valid triangular 
meshes. Hence, for any triangle t with edge lengths a^^\ h^^\c^^^ in p(^) and a^^\ h^^\c^^^ in P^^\ the 
following inequalities hold: 

a(0) + 6(0) > c(0), 6(0) + c(0) > a(0), c^o) + a^o) > 6(o) 

and 

ad) + 6(1) > c(i), Ml) + c(i) > a(i), c(i) + a^^^ > b^^l 

In P(*), a(*) = (1 - i)a^°^ + ia^^^, M*) = (1 - t)^^) + t6(i), c(*) = (1 - t)c(°) + fc(i) due to the linear 
interpolation of the end positions. Hence, a(*) + &(*) = (1 - i)a(°) + ta(i) + (1 - t)U-^'> + i6(i) = 
(l-t)(a(°) + 6(°)) + i(o(i) + 6(i)) > (l-i)c(o)+te(i). Similarly, we can show that 6(*)+c(*) > a(*) and 
c(*) + a(*) > &(*). Hence, p(*) is a valid triangular mesh. 

Part 2: The edge lengths of P(*) are Unear interpolations between the edge lengths of P^^) and p(i). 
Hence, the claim follows. 
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Part 3: When varying t continuously, the point p^^^ G 5 varies continuously. Hence, the coordinate 
of the lengths of all the edges vary continuously. Because a direction [sin v cos sin v sin u, cos v]^ varies 
continuously if u and v vary continuously, the normal directions vary continuously. Because all the vertex 
positions of the mesh p(^) are uniquely determined by continuous functions of those quantities, all vertex 
positions of P^^^ vary continuously. □ 

Note that we do not need to solve minimization problems to find the shortest path in shape space as 
in Kilian et al. ||9l. The only computation required to find an intermediate deformation pose is a graph 
traversal of L)(M). 

Because D{M) has complexity 0(n), we can traverse D{M) in 0(n) time. Hence, we can compute 
intermediate deformation poses in 0(n) time each. We denote this the polygon algorithm in the following. 



4 Theory of Shape Space for Skeleton Morphing 

A very similar shape space to the one presented in Section |3] can be used to isometrically morph between 
two topologically equivalent skeletons. Let a skeleton in be a set of joints connected by links arranged 
in a tree-structure. That is, we can consider a skeleton to be a tree in consisting of n vertices and n — 1 
edges. 

The shape space presented in Section [s] can be simplified to a shape space for skeletons in R^ with the 
property that interpolating linearly in shape space corresponds to the most isometric morph in R^. The 
dimensionality of the shape space is linear in the number of links of the skeleton. 

We start with two skeletons iS^^^ and 5^^^ corresponding to two near-isometric poses. We assume that 
the point-to-point correspondence of S^^"^ and 5^^) are known. Hence, we know the tree structure T with 
two sets of ordered vertex coordinates y^^^ and V^^"^ in R^. As before, we first find the best rigid alignment 
of the two skeletons. 

The skeleton shape space Ss is defined in a similar way as S. We assign an arbitrary but fixed root to 
T and traverse the edges of T in a depth-first order. We assign an arbitrary order to the edges incident on 
each vertex of T. The first 3 coordinates of s correspond to the coordinates of the root in T. The next n — 1 
coordinates of s are the lengths of the edges in T in depth-first order. The final 2(n — 1) coordinates of s 
describe the unit directions of the edges in spherical coordinates in depth-first order. All edges are oriented 
such that they point away from the root. Note that the shape space S has dimension 0(n). 

Interpolating linearly between points in Ss is performed the same way as interpolating linearly between 
points in S. Namely, edge lengths are interpolated linearly and unit directions are interpolated via SLERP. 
With a technique similar to the proof of Theorem [T] we can prove the following theorem. In the proof, we 
do not need to consider a dual graph, but we can simply traverse the tree T in depth-first order to propagate 
the information. 

Theorem 2 Let T be the underlying tree structure of the skeletons S^^^ and S^^\ The linear interpolation 
5^^) between s^^^ and s^^^ in shape space Ssfor < t < 1 has the following properties: 

1. The skeleton S^'^^ G R^ that corresponds to s^'^^ G Ss is uniquely defined and has the underlying tree 
structure T. We can compute this tree using a depth-first traversal of the tree in 0{n) time. 

2. IfS^^^ and S^^^ are isometric, then S^^^ is isometric to andS^^\ IfS^^^ andS^^^ are not perfectly 
isometric, then each edge length ofS^^^ linearly interpolates between the corresponding edge lengths 
ofS^^^ and S^^\ 

3. The coordinates of the vertices ofS^^^ are a continuous function oft. 
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This theorem allows us to morph isometrically between the skeletons of two shapes corresponding to 
two postures of the same articulated object in 0{n) time. Figure [2] shows an example of such a morph. 




(a) (b) (c) 

Figure 2: Example of isometric morph between two skeleton poses. Input poses are shown in (a) and (c) 
and morph for t = 0.5 is shown in (b). 

In the remainder of this paper, we will focus our attention on morphing between triangular meshes. 

5 Generalization to Triangular Meshes 

This section extends the shape space S from Section [3] to arbitrary connected triangular meshes. However, 
we can no longer guarantee the properties of Theorem [T] because the dual graph of the triangular mesh M 
is no longer a tree. 

Given two triangular meshes S^^"^ and 5^^) corresponding to two near-isometric poses of the same non- 
rigid object with known point-to-point correspondence, we know one mesh structure M with two sets of 
ordered vertex coordinates V^^"^ and V^^"^ in R^. 

As before, we can use this information to find the best rigid alignment in R^. We do this before 
representing the shapes in a shape space. As outlined above, this alignment has a major influence on the 
result of the morph. 

As before, we can represent S'(^) and S^^"^ as points 5^^^ and 5^^^ in a shape space S using the same 
shape space points as in Section [s] Let 5^^^ be the linear interpolation of 5^^^ and 5^^^ in 5, where the 
linear interpolation is computed as outlined in Section pi The existence of a mesh S^^"^ G R^ that has the 
underlying mesh structure M and that corresponds to is no longer guaranteed. This can be seen using 
the example shown in Figure [s] Figure [sja) and (b) show two isometric meshes S^^"^ and S^^\ The dual 
graph D{M) of the mesh structure M is a simple cycle. Note that although the start and the end pose 
are isometric, we cannot find an intermediate pose that satisfies all of the interpolated normal vectors with 
SLERP and that is isometric to S^^^ and S^^\ 




(a) (b) 

Figure 3: Example of isometric triangular meshes where intermediate poses interpolating all normals and 
edge lengths do not exist. 

Let M consist of n vertices. As M is a planar graph, M has 0(n) edges and 0(n) triangles. The shape 
space S is defined using the same shape space points as in Section [3] The shape space S has dimension 
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0(n). As before, we interpolate linearly in shape space by interpolating the edge lengths by a simple linear 
interpolation. 

Observation 1 Given a triangular mesh S^^"^ with underlying mesh structure M, point 5^^^ in S is uniquely 
determined. However, the inverse operation, that is computing a triangular mesh S^'^^ given a point s^'^^ G 
S, is ill-defined. 

This is illustrated in Figure |3] 

To compute a unique triangular mesh S^^^ given a point s^^^ G S that linearly interpolates between s^^^ 
and s^^\ such that S^'^^ approximates the information given in 5^^^ well, we use the dual graph D{M) of 
M. Unlike in Section [3} D{M) is not necessarily a tree. Our algorithm therefore operates on a minimum 
spanning tree T{M) of D{M). The tree T{M) is computed by assigning a weight to each arc e of D{M). 
The weight of e is equal to the difference in dihedral angle of the supporting planes of the two triangles 
of M corresponding to the two endpoints of e. That is, we compute the dihedral angle between the two 
supporting planes of the two triangles of M corresponding to the two endpoints of e for the start pose 5^^) 
and for the end pose S^^\ respectively. The weight of e is then set as the difference between those two 
dihedral angles, which corresponds to the change in dihedral angle during the deformation. The weight can 
therefore be seen as a measure of rigidity. The smaller the weight, the smaller the change in dihedral angle 
between the two triangles during the deformation, and the more rigidly the two triangles move with respect 
to each other. As T{M) is a minimum spanning tree, T{M) contains the arcs corresponding to the most 
rigid components of M. 

Similar to Section [s} we compute S^^^ by traversing T{M). However, unlike in Section [s} setting the 
vertex coordinates of a vertex v of S^'^^ using two paths from the root of T{M) to two triangles containing v 
can yield two different resulting coordinates for v. An example of this situation is given in Figure |4} where 
the coordinates of v can be set by traversing the arcs 62, and 63 of T{M) or by traversing the arcs ei, 64, 
and 65 of M starting at root{T{M)). We call the different coordinates computed for v in T{M) candidate 
coordinates of v. Our algorithm computes the coordinates of each vertex v G S^^^ as the average of all the 
candidate coordinates of v. 



V 




Figure 4: A mesh M with its dual minimum spanning tree T{M). 

Let us analyze the maximum number of candidate coordinates that can occur for a vertex in S^^\ Let e 
denote an edge of M such that D{e) is in T(M). Let v denote the vertex of 5^^^ opposite e in the triangle 
corresponding to an endpoint of D(e), such that the coordinates of v are computed when traversing D(e). 
This situation is illustrated in Figure [5] Let di and d2 denote the total number of candidate coordinates of 
the two endpoints of e. By traversing D{e), we compute did2 candidate coordinates for v. We can therefore 
bound the number of candidate coordinates of v computed using the path through D{e) by did2. Note that 
the number of candidate coordinates for the two endpoints of the first edge is one. Furthermore, each vertex 
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V can be reached by at most deg{v) paths in T{M), where deg{v) denotes the degree of vertex v in M. As 
each path in T{M) has length at most m — 1, where m = 0(n) is the number of triangles of M, we can 
bound the total number of candidate coordinates in S'(^) by J2vev '^^~^deg{v) = 2n2^~^, where V is the 
vertex set of M. 




Figure 5: Illustration of how to bound the number of candidate coordinates ofv computed using the path 
through D{e). 

Our algorithm finds a triangular mesh S^^^ corresponding to 5^*) that is isometric to S^and^W if such 
a mesh exists, because all of the candidate coordinates are equal in this case and taking their average yields 
the desired result. If there is no isometric mesh corresponding to s^'^\ our algorithm finds a unique mesh 
that weighs all the evidence given by T{M) equally. By choosing T{M) as a minimum spanning tree based 
on weights representing rigidity, we allocate rigid parts of the model more emphasis than non-rigid parts. 
The reason for this is that in most near-isometric morphs, triangles close to non-rigid joints are deformed 
more than triangles in mainly rigid parts of the model. We conclude with the following Lemma. 

Proposition 1 Let S^^^ and S^^^ denote two isometric connected triangular meshes and let s^^^ and s^^^ 
denote the corresponding shape space points, respectively. We can compute a unique triangular mesh 
S^'^^ approximating the information given in the linear interpolation s^^^O < t < 1 of s^^^ and s^^\ in 
exponential time. We find a triangular mesh S^'^^ corresponding to s^'^^ that is isometric to 5(0) and 5(1) if 
such a mesh exists. 

The algorithm can easily be extended to work for a non-connected triangular mesh M by removing rigid 
transformations for each connected component of M using local coordinate systems. We can then adapt 
the algorithm by finding the dual graph D{M) and a minimum spanning tree T{M) for each connected 
component of M. With this information, we can traverse the graph as described above. We denote this the 
exponential algorithm in the following. 

6 Efficient Algorithm to Deform Triangular Meshes 

As the exponential algorithm is limited for pragmatic reasons to triangular meshes with few vertices, this 
section describes a more computationally efficient algorithm to find the deformed poses. 

To compute a unique triangular mesh S'(^) given a point 5(^) G S that linearly interpolates between 5(0) 
and 5(1), such that S'(^) approximates the information given in 5(^) well, we use the minimum spanning tree 
T{M) of the dual graph D{M) of M as before. An approach that reduces the total number of candidate 
coordinates of each vertex of S'(^) to 0{n) is to restrict to one the number of times each edge of T{M) can 
be traversed. We traverse T{M) in depth-first order. When an edge D{e) is traversed, we add candidate 
coordinates to one vertex v as shown in Figure |5] However, we only add at most a linear number of 
candidate coordinates to v as described below. Denote the vertices of e by vo{e) and vi{e){e), denote the 
number of candidate coordinates that were added to V(){e) and vi{e), respectively, during the traversal of 
T{M) before traversing the edge D{e) by di and 0^2, and let VQ{e) be the vertex of e that was updated more 
recently in the traversal of T{M). Let the candidate coordinates of V{){e) {vi{e), respectively) be given by 
c\, . . . ,cf^ (c2, . . . , C2^, respectively) ordered from the least recently to the most recently added candidate 
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coordinate. When traversing D{e), we add d2 candidate coordinates to v by computing coordinates of v 
based on the candidate pairs {cf^ , C2), . . . , {cf^ ^c^^)- 

This strategy computes 0{n) candidates per vertex of S^^\ and hence a total of O(n^) candidates, 
thereby avoiding the computation of an exponential number of candidate coordinates. To find the final 
coordinate of a vertex v, we average all of the candidate coordinates of v. We denote this the averaging 
algorithm in the following. 

The averaging approach was found to yield satisfactory results in all test cases. Just as in the algorithm 
by Kilian et al. (91, the results depend heavily on the initial rigid alignment of the shapes. 

Furthermore, the results of the averaging algorithm depend on the order of the vertices and edges in M. 
Recall that the first vertex and its first incident edge are known. As we limit the number of times an edge of 
T(M) can be traversed, not all of the information is propagated through the graph. Hence, the first vertex 
and the first edge influence the final result. We can eliminate this dependence at the cost of a higher running 
time by taking 0{n) shape spaces; one with each possible oriented edge as first edge. The running time of 
this algorithm is O(n^), which is too high for practical applications. We therefore can also use a smaller 
number of first edges that are found using Voronoi sampling [7|. If 0(log n) samples are used, the running 
time of the algorithm becomes 0{n log n) and the influence of the first edge is expected to be low. We do 
not include experiments obtained using these extensions of the approach as the averaging algorithm yields 
similar results as the less efficient algorithms in all of our experiments. 

7 Experiments 

This section presents experiments using the algorithms presented in this paper. We show exact morphs of 
triangulated 3D polygons, morphs of models with few vertices computed using the exponential algorithm, 
and morphs of larger triangular meshes using the averaging algorithm. 

The experiments were conducted using an implementation in C++ on an Intel (R) Pentium (R) D with 
3.5 GB of RAM. OpenMP was used to improve the efficiency of the algorithms. To compute the minimum 
spanning tree T(M), the boost graph library |[T6l was used. 

7.1 Deforming Triangulated 3D Polygons 

We demonstrate the efficient polygon algorithm using the simple polygon shown in Figure]?] We deform the 
polygon shown in Figure [7] (a) to the polygon shown in Figure |7](i). The morph is illustrated in Figures [7] 
(b)-(h). All of the intermediate poses are isometric to the start and end poses. The overlayed poses are 
shown in Figure [6j The running time of this example is less than 1 second. 




Figure 6: Isometric morph of a simple polygon. The start polygon is a 3D polygon obtained by discretizing 
the curve y — sin(x) and by adding thickness to the curve along the z-direction. The end polygon is 
similarly obtained from y — — sin x. 

7.2 Deforming General Triangular Meshes 

We present experimental results for the exponential algorithm and the averaging algorithm. First, we run 
the exponential algorithm on one model with few vertices to demonstrate the quality of the results. The 
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(a) (b) (c) (d) (e) (f) (g) (h) (i) 



Figure 7: Isometric morph of a simple polygon from pose (a) to pose (i). Intermediate poses obtained using 
the polygon algorithm are given from left to right. 



model we use to test the approach is shown in Figure [8j We aim to smoothly and isometrically deform 
the pose shown in Figure [8|a) to the pose shown in Figure [8|i). As mentioned above, there is no isometric 
deformation between the poses that interpolates the triangle normals. The result of our algorithm is shown in 
Figures[8jb)-(h). Note that all triangle normals are interpolated and the symmetry of the model is preserved. 
Furthermore, all edge lengths with the exception of the edges of the four top faces are interpolated. 

(a) (b) (c) (d) (e) (f) (g) (h) (i) 

Figure 8: Isometric morph of a cycle from pose (a) to pose (i). Intermediate poses obtained using the 
exponential algorithm are given from left to right. 



Second, we demonstrate the quality and efficiency of the averaging algorithm. The first experiment 
morphs between animated poses of the armadillo model. The models are chosen from the AIM ©SHAPE 
repository Q The models contain 331904 triangles and 165954 vertices. Testing our algorithms on these 
models emphasizes the practical use of our method. The results are shown in Figure |9] Note that the 
intermediate poses are visually pleasing. Although some small bumps appear on the right arm of the 
armadillo in Figure [9] (b) and on the left arm of the armadillo in Figure [9] (e), most details are preserved 
during the morph. 

Figure [To] shows the morph between two poses of a model of a human being. The deformed pose of 
the model was found using the automatic technique by Baran and Popovic [3|. The models contain 10002 
vertices. The two given meshes are shown in Figures [TOl^a) and (i). Intermediate poses using the averaging 
algorithm are shown in Figures [TOjb) to (h). Note that although small bumps appear on the leg of the 
morphed surface, the morphs are intuitive. 

Furthermore, we morph between heads from the CAESAR database 1 13 1. The correspondence between 
the two head models was found using the approach by Xi et al. [22]. The models contain 11102 vertices. 
The two given meshes are shown in Figures [TTJa) and (i). Intermediate poses obtained using the averaging 



algorithm are shown in Figure 11 (b) to (h). Note that the morphs are visually pleasing. 



Finally, we morph between a few poses of the Alien model shown in Figure [T2j The Alien model 
is chosen from the Princeton Shape Benchmarlsj^ and animated to obtain multiple postures with known 
correspondences using the automatic technique by Baran and Popovic O. All of the Alien models contain 
429 vertices. The results are shown in Figure[T2]and we can see that all of the morphs are visually pleasing. 

For all of the experiments conducted, we measured the time efficiency to compute one intermediate 
pose at t = 0.5 and the energy 



Q = \^{\\vo{e) - vr{e)\\ - l{e)y 

eeE 



^ http : //shapes . aimatshape .net/releases .php 
^http ://shape.c s .princeton.edu/benchmark/ 
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(d) (e) (f) 

Figure 9: Isometric morph of the armadillo model from pose (a) to pose (c) and from pose (d) to pose (f). 
The intermediate poses obtained using the averaging algorithm for t = 0.5 are given in figures (b) and (e). 

where E is the set of edges of M, vo{e) and vi{e) are the vertices of e in S^^\ ||a, 6|| denotes the Euclidean 
distance between a and b, and /(e) is the desired length for e stored in s^^\ The running times and the 
quality of the results are summarized in Table [T| Note that the developed implementation is non-optimized 
and experimental. The running times might be improved by implementing some of the algorithms on the 
GPU. 

8 Conclusion 

We presented an approach to morph efficiently between isometric poses of triangular meshes in a novel 
shape space. The main advantage of this morphing method is that the most isometric morph is always found 
in linear time when triangulated 3D polygons are considered. For general triangular meshes, the approach 





Alien Row 1 


Alien Row 2 


Human 


Head 


Armadillo Row 1 


Armadillo Row 2 


n 


429 


429 


10002 


11102 


165954 


165954 


time (sec) 


< 1 


< 1 


8 


10 


1813 


1763 


Q 


0.000719 


0.000064 


0.042674 


0.003583 


12277.92 


70637.80 



Table 1: Summary of the experimental results. 
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(a) (b) (c) (d) (e) (f) (g) (h) (i) 

Figure 10: Isometric morph of a human from pose (a) to pose (i). Intermediate poses obtained using the 
averaging algorithm are given from left to right. 




(a) (b) (c) (d) (e) (f) (g) (h) (i) 

Figure 11: Isometric morph of a head from pose (a) to pose (i). Intermediate poses obtained using the 
averaging algorithm are given from left to right. 

cannot be proven to find the optimal solution. However, this paper presents an efficient heuristic approach 
to find a morph for general triangular meshes that does not depend on solving a non-linear optimization 
problem. 

The presented experimental results demonstrate that the heuristic approach yields visually pleasing 
results. The approach is not invariant with respect to the order of the vertices of the mesh, but can be 
modified to have this property at the cost of a higher running time. 

An interesting direction for future work is to find an efficient way of morphing triangular meshes while 
guaranteeing that no self-intersections occur. For polygons in two dimensions, this problem was solved 
using an approach based on energy minimization |[8ll . 
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